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Abstract 
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which allow for explicit breather solutions. Especially we obtain exponentially 

localized solutions for <i-dimensional lattices with d = 2,3. 
02.10.Gd, 05.45.Xt, 05.45.Yv, 63.20.Pw, 63.20.Ry 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



The understanding of dynamical localization in classical spatially extended and ordered 
systems |l[ experienced recent considerable progress @, ||, |4|]. Specifically time-periodic 
and spatially localized solutions of the classical equations of motion exist, which are called 
(discrete) breathers, or intrinsic localized modes. The attribute discrete stands for the 
spatial discreteness of the system, i.e. instead of field equations one typically considers the 
dynamics of degrees of freedom ordered on a spatial lattice. The lattice Hamiltonians are 
invariant under discrete translations in space. The discreteness of the system produces a 
cutoff in the wavelength of extended states, and thus yields a finite upper bound on the 
spectrum of eigenfrequencies fl q (phonon band) of small-amplitude plane waves (one usually 
assumes that usually for small amplitudes the Hamiltonian is in leading order a quadratic 
form of the degrees of freedom). If now the equations of motion contain nonlinear terms, 
the nonlinearity will in general allow to tune frequencies of periodic orbits outside of the 
phonon band, and if all multiples of a given frequency are outside the phonon band too, 
there seems to be no further barrier preventing spatial localization (for reviews see 0, ||). 

Discrete breathers have been recently experimentally detected in weakly coupled waveg- 
uides 0, MX solids P and Josephson junction ladders |7|]. The broad spectrum of applica- 
bility of the localization concept described above makes it worthwhile to continue theoretical 
efforts to characterize the properties of discrete breathers. 

So far there is very little knowledge about explicit forms of breather solutions. Except 
for trivial limiting cases like the antiintegrable limit, i.e. the case of uncoupled degrees of 
freedom ||, we know only about the solutions of the Ablowitz-Ladik lattice || (an inte- 
grable one- dimensional variant of the nonlinear Schrodinger equation). What is generically 
available is the abstract knowledge about existence or nonexistence of discrete breathers for 
a specific system, and the spatial decay properties far from the center of the breather, where 
due to the smallness of the amplitudes linearizations or other perturbation techniques are 
applicable. Note that generically breathers can appear in nonintegrable systems. 
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From the above it appears that nonintegrability spoils the possibility to obtain analytic 
forms of the breather solution. We will show that this is not the case by constructing Hamil- 
tonians which allow for explicit solutions and are most probably not integrable. Moreover we 
will even construct solutions for two- and three-dimensional lattices. Although these models 
are not motivated by certain applications, the study of their properties can be helpful with 
respect to discrete breathers. 



II. A BOND-ORDERED QUASI-LINEAR CHAIN 

In this section we consider a one-dimensional model with the Hamiltonian 

H = w h + w p , w k = Y,\p 2 i . % = ^£(z/-W^) , *({*!'}) = — + — • 

{ l l x Xi + i Xi 

(1) 

The integer I marks the lattice sites of the chain. The equations of motion read 

xi=Pi , Pi = — ^— • (2) 

OXi 

Since we are interested in obtaining solutions which decay to zero at spatial infinities and 
can be interpreted as excitations above some regular ground state x\ = 0, we demand that 
h(s) behaves regularly and especially that h(s — > ±oo) does not diverge. The potential 
energy W p is a homogeneous function of the coordinates xi since W p ({Xxi}) = X 2 W p ({xi}) . 
The homogeneity of the potential function can be used to separate time (G(t)) and space 
(ui) variables as done e.g. in |H|, [ |TT|j : 

x l (t)=u l G(t) , G=-kG , -K Ul = - dWp tt Ul ' }) . (3) 

Here k > is needed to ensure boundness of the solutions. Indeed the time dependence is 
then given by 

G(t) = Acos(ujt + 4>) , uo 2 = k . (4) 
The equations for the spatial amplitudes ui read 



K=(l-^±i) h ( Sl ) + (l-^\>( Sl ) Pl 

( 1 - —) h(s^) -(l-—Y tffa-OPw , (5) 



Ui J \ Ui 



Ul u l+1 
Pi = 6 

Ul+l Ui 



and si = si({ui'}). Here f'(x) means the first derivative of / w.r.t. x. We are looking for a 
solution of {ui>} which is localized in space, i.e. ui^± OQ — > 0. In order to find such a solution 
to (|5]) we assume that si is constant in the spatial tails. This condition is equivalent to 
having exponential decay. Taking 

u, = (-l)'e-^l (7) 

and combining the two different cases I = (center of the solutions) and I ^ (tails of the 
solution) we find 

h(-s) = -s{s + l)h'{-s) , s = 2cosh(/5) (8) 

and 

' 1+ 7rr( 1+ f (1 -^) 



k= (2 + s(l - 7)) h{-s) 



- (-V • (9) 



Suppose there exists a value for s = so > 2 such that (|8]) is satisfied. If for this value we 
also have h(—s) > then @ defines a positive value for k. Moreover the found solution so 
of (H) will be structurally stable against changes in h(s). Such models can be generated by 
starting with the trivial case h(s) = 1. A strong enough local distortion in h(s) at s < —2 will 
generate solutions of the above equations, preserving the overall positivity and boundness 
a > h(s) > for all s with finite a. Thus the state x\ = c (here c is an arbitrary constant) will 
be a state with minimum energy H = 0, with all other trajectories having larger energies. We 
can therefore interprete the found explicit localized time-periodic solution (|7|) as a discrete 
breather solution, an excitation above a classical homogeneous ground state. 
Let us consider an example. Choose 

h(s) = l + ae~ bs2 . (10) 



For parameters 

e 3/2 3 

a> , 6>- (11) 

2 + 3^3/(26) § 

there is at least one solution to ([5]) with s > 2. This solution is structurally stable against 
perturbations in h(s). E.g. for b = 0.5 and a = 0.7 the solution is s ~ 2.06781, /3 ~ 0.25967. 

Why did we choose a staggered solution u\ ~ (— l) 1 ? This is motivated by the fact that 
we consider perturbed harmonic chains (h(s) = 1) for which the spectrum of plane waves is 
acoustic. Thus discrete breathers in order to localize should show up with frequencies above 
the acoustic phonon band which implies out-of-phase motion of nearest neighbours, or simply 
staggered solutions. It is not that simple to apply these arguments to the case of h(s) ^ 1, 
since there is no simple way to linearize the equations of motion around xi = in the general 
case. Without going into further details let us state here, that using a nonstaggered ansatz 
one arrives at equations which do not ensure positivity of k in general. 

Let us study the problem of small amplitude excitations above the groundstate of ([[]). 
First we recall that we consider only positive and bounded functions h(s). Then there exists 
a continuous family of ground states, i.e. solutions with lowest possible energy E — and 
x\ = which are given by 

xi = c . (12) 

Notice that the Hamiltonian ([!]) is not invariant under transformations x\ — > c + xi, yet the 
ground state energy is degenerate. An expansion around one of the ground states yields 

W P = % E [M2) (5/ - <W 2 + h'(2) (St - S l+1 ) A + 0(S 5 )} , S t = ^ - 1 • (13) 
1 i ' ' 



c 



The neglected terms of fifth and higher orders in ([D^) are not invariant under transformations 
Si —>■ c + Si. However the terms up to fourth order are invariant. Taking into account only 
terms up to fourth order thus yields a so-called Fermi-Pasta-Ulam chain for the dynamics 
of small deviations from the ground state x\ = c. Note that we can not simply perform the 
limit c — > since the expansion (|I3|) is valid only if |^| c. The above found breather 



solution, which decays to zero at infinities (and not to c 7^ 0), can not be easily deformed in 
order to decay to c 7^ at infinities. However it is wellknown that (|i3|) allows for discrete 



breather solutions if h'{2) > (which can not be given in a closed analytical form) (e.g. |T2 



|T3|j , |H|). So we conclude that for the considered model (|1|) discrete breather excitations 
above the ground state x\ = c 7^ are solutions with generic features, and the ground state 
x\ = allows for discrete breather excitations given in a closed analytical form. 

In the last part of this section we will discuss the existence of different variants of discrete 
breathers with £;_^oo — > asymptotics. The existence of the above derived discrete breather 
solution can be interpreted as follows. Our ansatz ui = (— l)^ - ^'^ contains together with 
the parameter k (coming from separating time and space in ([3D) two parameters to be 
determined - (3 and k. With our ansatz we found two equations - one for the spatial wings 
of the solution, and one for the center. Two equations with two variables can be solved in 
general, and the additional inequality k > will serve as an additional restriction for the 
choice of possible functions h(s), but will not change the fact that once solutions are found, 
they will be in general structurally stable against changes in h(s). Let us now look for a 
solution of the form u\ = e~@ l for / > 0, u\ = ae^' 1 for I < — 1. We now have four equations 
- two in the tails, one at Z = and one at / = —1, and four parameters - K,,j3,j3',a. So 
we conclude that in general such solutions will exist. Indeed, the solution from above is a 
variant of the more general case discussed here with a = 1 and f3' = j3. So we can expect 
in general a countable set of other solutions with a ^ 1 and f3' 7^ j3. All these solutions 
will have a closed analytical form with parameters to be determined numerically from the 
mentioned equations. 

III. SITE ORDERED MODELS IN D DIMENSIONS 

Consider a d- dimensional hypercube lattice with a scalar coordinate x\ associated to 
each lattice site. The site index 1 = ■■■,ld) is a rf-dimensional vector with integer 

components /j. Consider the operator L defined as 
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Lxi = x v ■ 

i',|i'-i|=i 



(14) 



Here |1| 2 = l( + l* + ... + lj. Also define 



si 



(15) 



The Hamiltonian is given by a sum over kinetic and potential energies as in ([[]), with the 
potential energy 



(16) 



Assuming time-space separability as in (§) we obtain 



G = -kG 



(17) 



for the time dependence. Again we need k > to ensure bounded solutions. 

The spatial coordinates have to satisfy an equation similar to (|3|). Let us calculate the 
derivative 



dW p d 



If \ 2 

2 h ( s h 



du\ du\ 

Evaluation of these equations leads to the result 

1 



— L 



2 ( Lx iJ K s h 



(18) 



KU\ = —-Uisfh'(si) + L [uiSih(si)] + —L u\s{h\s\) 



(19) 



A closer study of equation (]19f) shows that it could be easily solved if s\ is essentially constant 
everywhere on the lattice, more precisely everywhere except for one site 1 = 0: 



Si^o — s l j s i=o — s o 



In that case evaluation of flT9"p yields 



K = slh(si) + 



Mo 
Ml 



s h(s ) + 2 S o h '( s o) 



k = slh(si) , |1| > 1 

i .... .i 



k = s /i(so) , 1 = 



(20) 



(21) 
(22) 
(23) 



Defining a new function 

g(a) = \s 2 h{s) (24) 

equations (pT])-(p3|) are reduced to 

K = 2g(s ) , flf(s ) = flf(si) , #'(s ) = fl^(si) • (25) 



These relations (|25[) are in fact conditions on the choice of Hamiltonians, i.e. of the function 
h(s). They are so far rather general, but we have to specify the values sq, s\. These values 
will be connected with the solution of (|i~9D through the conditions (|20|) , which actually 
constitute a linear equation: 

Lu\ = S1M1 . (26) 

A solution to this equation can be cast into the form 

f, , , cos(Zigi)cos(Z 2 g2)-cos(Z d g d ) 

u x = / d qi dq 2 ..Aq d — — — — r— — , 27 

J 2 (cos(gi) + cos(g 2 ) + ••• cos(g d )) - s 1 

where the integration extends for each variable from — it to tt. Then the value for so is 
given by 

' f d qi dq 2 ...dq d — -— -— . (28) 



s -st J "' 2 (cos(gi) + cos(g 2 ) + ••• cos(q d )) - s x 

In other words, given the solution fl27|) to ©, we can generate the corresponding HamiL 



tonian by solving (|25D with the additional constraint (^3), which fixes the value sq relative 
to si. The function g(s) can then be constructed in the following way: first choose Si, 
then determine sq, finally find a function g(s) which is positive in si and whose values and 
derivatives are equal in both so and s^. This function g(s) will then define h(s) with (^ 
and thus the Hamiltonian with potential energy ([16]). 

Notice that (|27|) is exponentially localized around 1 = if |si| > 2d. 
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IV. CONCLUSIONS 



We obtained explicit discrete breather solutions for several classes of Hamiltonian lattice 
models with different lattice dimensions. The outlined approach can be extended to other 
classes of lattice models with homogeneous potentials. 

An interesting question for future studies will be the linear stability analysis of the 
obtained solutions, which is also closely related to the spectrum and character of small 
amplitude fluctuations around the state x\ = 0. Note that any solution which satisfies time- 
space separability can be continuously tuned to X\ = by letting A — > in (||), regardless 
of its spatial profile. 

Finally one can extend the solutions of the site-ordered models to multisite breathers, 
for which 

si^i m = , si m = s m (29) 

where the number of 'defect' sites l m is finite and the largest distance between any pair of 
'defect' sites is finite. Work is in progress. 

1 Permanent address: Institute of Chemical Physics, Kosygin St. 4, 117334 Moscow, Russia. 
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